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Abstract. Aspects of turbulence in protostellar accretion discs are being reviewed. 
. The emergence of dead zones due to poor ionization and alternatives to the magneto- 

rotational instability are discussed. The coupling between dust and gas in protostellar 
accretion discs is explained and turbulent drag is compared with laminar drag in the 
Stokes and Epstein regimes. Finally, the significance of magnetic ŕield generation in 
turbulent discs is emphasized in connection with driving outfiows and with star-disc 
coupling. 
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í> ' 1. Introduction 

O 

. 

. One of the amazing properties of accretion disc theory is the fact that it is applicable 
over a vast range of scales, from discs around supermassive black holes to those 
around stellar mass black holes and around neutrón stars, to young stellar objects 
and possibly even protoplanets. This is mainly because the equations of hydrodynamics 
and magnetohydrodynamics are scale-invariant over broad ranges. In all cases accretion 



is caused by the action of the Reynolds and Maxwell stresses. The stress is assumed 
to be expressible in terms of the radial angular velocity gradient multiplied by some 
turbulent viscosity coefficient. However, a major problém is that in protostellar discs 
the temperatures are rather low, so the degree of ionization depends significantly on 
cosmic ray ionization. At a radial dištance of about 1 AU the regions near the midplane 
are strongly shielded, which implies the possibility of "dead zones" (Gammie 1996) 
where magnetic effects cannot play a role. This leads to our first topic. 

Accretion is necessary for accumulating the mass of the centrál object which is 
in this case a protostar. This accretion is presumably of turbulent origin, because the 
molecular viscosity is obviously far too small. Nevertheless, the molecular or microscopic 
viscosity together with microscopic magnetic diffusion are crucial, and they are in fact 
the main processes that can heat the disc-in addition to irradiation from the centrál 
star onto the our surface of the disc. However, it is still unclear whether in a turbulent 
disc planét formation occurs because of turbulence or in spite of it. This will be our 
second topic. 
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Finally, we address the issue of outflows from the star and from the disc. The 
occurrence of collimated outflows in star-forming regions has always been associated 
with magnetic fields. Given that magnetic fields can also be produced within the disc 
from which such outflows emanate, the question emerges whether the field necessary for 
collimation has to be an external interstellar field, or whether the disc field suflices. In 
addition, there will be the stellar field and its coupling to the ambient magnetic disc. 
This is thought to be critical for spinning up the protostar during its formation and 
spinning it down later during its evolution. This is the final topic discussed in this 
review. We štart however with a discussion of the mean accretion stress, its definition, 
and how it is used. 



2. The mean accretion stress 



In order to get an idea about the náture of turbulence in accretion discs, it is useful 
to consider simulations in a local cartesian geometry with shearing sheet boundary 
conditions in the radial direction. An important output parameter of such simulations 
is the Shakura-Sunyaev viscosity parameter, a S s> which is a non-dimensional measure 
of the turbulent viscosity u t , in terms of the sound speed c s and the disc scale height H, 
i. e. 

ľt = ®ssc s H. (1) 
This coeflicient is used analogously to the molecular viscosity and hence is used to 
multiply the velocity gradient to get the accretion stress. The main difference is however 
that the turbulent viscosity multiplies the gradient of the mean velocity while the actual 
viscosity multiplies the actual turbulent velocity gradients, which can be much higher 
locally on very small scales. In the mean field description the mean stress is written as 

dU 

Tzucf, = -pV t VO — , (2) 

where w is the cylindrical rádius, p is the density, fž is the mean angular velocity, 
and ~T m ^ is the total 'horizontál' stress. This component of the stress tensor enters the 
conservation equation for the mean angular momentum, 

^ (pzu 2 U) + V • (pUU^ + «7r w ,) = 0. (3) 

Here, vj = (1, 0, 0) T is the unit vector in the direction away from the axis, and overbars 
denote mean quantities that are obtained by averaging over the azimuthal direction. 
The velocity is then written as U = U + u, and fluctuations of the density p are 
ignored. Also, for simplicity, we have ignored the vertical component of the stress, t z ^. 
Roughly, the two terms under the divergence balance, so a positive stress, r ro< ^ > 0, 
causes a negative (inward) accretion velocity, < 0. 

In a theory that works with mean velocities, U, as opposed to the actual one, U, 
one needs to calculate the contributions of small scale velocities and magnetic fields on 
the mean flow through the sum of mean Reynolds and Maxwell stresses, 



-b zu b < p/p + pu^Uý, (4) 



Turbulent protostellar discs 



3 



0.020 



0.015 



ä o.oio 



0.000 

















: 

, i 




_ 


! 












U. 

: : | 


_ 




fcAí 


|, 










'■'.'! : 








f *1 




■í 

.-i' „ 




















i/l'* 








-'„; A \ / 


í i \\ 












,/' 'K', 


/ — 

/ť * \ 'V ~h* 





















150 



200 
t/T rot 



250 



Figúre 1. Time serieš of the Shakura-Sunyaev viscosity alpha, ass(í) (dotted line), 
compared with a fit of the form equation ((6]) in tcrms of the instantaneous magnctic 
field strength (solid line). Time is normalized in terms of the orbital time, T rot = 2tt/H, 
where fž is the local angular velocity. Note that ass fiuctuates strongly in time about 
an average value of around 0.007. [Adapted from Brandenburg (1998).] 



where b and u denote the fluctuating components of the magnetic and velocity fields, 
and Hq is the vacuum permeability. The first term, i. e. the magnetic stress, is usually the 
largest term. The time dependence of the resulting stress determined from a simulation, 
and expressed in nondimensional form as 

ass(í) = ťW (pc s H\n) , (5) 
is shown in Fig. [TJ The simulation domain covers only the upper disc piane. The 
instantaneous time dependence follows closely that of the instantaneous magnetic field 
strength and can be described by a fit of the form 

f , s _ (0) . (B) (B) 2 . . 

where (B) is the mean magnetic field averaged over the entire upper disc piane, and 
Bq = ^fiopc s is the thermal equipartition field strength, with c s being the sound speed, 
and fit parameters af^ ~ 0.002 and atgg ~ 0.06 (Brandenburg 1998). Blackman et 
al. (2008) have shown that a similar relationship, albeit without the «gg term, can be 
established for a number of other simulations published in the literatúre. 

Simulations based on the magneto-rotational instability (Balbus & Hawley 1991, 
1998) show that the average value of ass(^) is small (around 0.007; see Hawley et al. 
1995, 1996, Stone et al. 1996, Brandenburg et al. 1995, 1996). Larger values of ass are 
obtained with an externally imposed vertical or axial magnetic field (e. g. Hawley et al. 
1995, Torkelsson et al. 1996), as is also predicted from equation d^J). However, a ss may 
well be smaller, depending on simulation detailsjij] Indeed, there is still an open issue 

í Here the scale height has been defined such that p ~ exp(— z 2 /H 2 ). Later we shall define a scale 
height such that p ~ exp(— \z 2 /H 2 ). Note that if H = H /y/2 were used in the defmition of ass 
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regarding the convergence of the numerical results (see, e. g., Fromang & Papaloizou 
2007, Fromang et al. 2007). 

The variation of «ss(^) shows a typical time scale of around 15 orbits. In a 
simulation with vertical density stratification this can be a consequence of the emergence 
of cyclic large scale dynamo action. The magnetic field reverses between two maxima, 
so the actual periód is around 30 orbits. Physically, this time scale is related to the 
turbulent diffusion time. This is now reasonably well understood and the details of 
the large-scale dynamo magnetic field depend obviously on boundary conditions. Such 
cycles are not expected in global geometry, as can be demonstrated by corresponding 
mean field calculations (Bardou et al. 2001). 

Global simulations of accretion disc turbulence have been studied by various groups 
since the late 1990s. A conceptionally simple approach is that of the cylindrical disc 
which omits vertical stratification. In such simulations the magnetic field is able to 
develop over large scales and can be substantially stronger than in local simulations, 
resulting therefore also in larger values of «ss; Armitage (1998) found ass ~ 0.1. Similar 
results where later confirmed with fully global simulations by Hawley (2000), who found 
values of ass in the range of 0.1-0.2. Simulations of global discs permeated by a magnetic 
field from a centrál star also obey equation (Q, but in the simulations of Steinacker & 
Papaloizou (2002) the magnetic field was weaker, resulting in ass values of around 0.004. 
Similar figures have also been reported by Nelson & Papaloizou (2003) and Fromang & 
Nelson (2006) in the context of turbulent protoplanetary discs with and without planet- 
disc interaction. Global simulations have now also been able to produce estimates for 
Reynolds and Maxwell stresses as a function of rádius using the disc-in-a-box approach 
(Lyra et al. 2008). Among other things they find that larger Mach numbers result in 
larger normalized accretion stresses, as quantified by the Shakura-Sunyaev viscosity 
alpha. 

3. Dead zones 

When the gas is cool, the main source of ionization is not thermal, but it is cosmic 
ray ionization from the galaxy and UV ionization from the centrál star (and possibly 
other nearby stars). Research in this direction has been pursued by Glassgold et al. 
(1997) and Igea et al. (1999). Cosmic ray ionization usually plays an important role 
at moderate and large distances from the centrál object and away from the midplane 
where cosmic rays are shielded. Here, the column density of the gas is so high that the 
degree of ionization is very small, and the magneto-rotational instability (MRI) cannot 
grow because of the large diffusivity of the field. In addition, magnetic fields and a 
magnetized wind from the centrál star have a tendency to shield the disc from cosmic 
rays. The resulting región of low ionization in the disc is therefore referred to as the 

in equation (TT]), the resulting value of ass would be larger by a factor %/2- Note also that in some 
definitions of ass the factor 3/2 from wdíl/dw — — | íž is neglected, in which case ass would be larger 
by yet another factor of 3/2, so altogether « 2.1. 
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dead zóne (Gammie 1996). 

Recent work of Fleming & Stone (2003) shows that, although the local Maxwell 
stress drops to negligible values in the dead zones, the Reynolds stress remains 
approximately independent of height and never drops below approximately 10% of the 
maximum Maxwell stress, provided the column density in that zóne is less that 10 
times the column density of the active layers. The non-dimensional ratio of stress to 
gas pressure is just the Shakura & Sunyaev (1973) viscosity parameter, a ss . Fleming 
& Stone (2003) find typical values of a few times 10~ 4 in the dead zones and a few 
times 10 -3 in the MRI-active layers. Similar results have recently been confirmed 
by Oishi et al. (2007). On the other hand, Inutsuka & Sano (2005) have questioned 
the very existence of dead zones, and proposed that the turbulent dissipation in the 
disc provides sumcient energy for the ionization. However, their calculation assumes a 
magnetic Reynolds number that is much smaller than expected from simulations (see 
Matsumura & Pudritz 2006). 

The radial extent of the dead zóne depends primarily on the disc column density 
as well as effects of grains and radiation chemistry (e. g., Sano et al. 2000, Matsumura 
& Pudritz 2005). Inside the dead zóne the magnetic Reynolds number tends to be 
below a certain critical value that is somewhere between 1 and 100 (Sano & Stone 
2002), making MRI-driven turbulence impossible. Estimates for the radial extent of the 
dead zóne range from 0.7-100 AU (Fromang et al. 2002), to 2-20 AU in calculations 
by Semenov et al. (2004). For smaller radii, thermal and UV ionization are mainly 
responsible for sustaining some degree of ionization (Glassgold et al. 1997; Igea et al. 
1999). The significance of the reduced value of ass in the dead zones is that it provides 
a mechanism for stopping the inward migration of Jupiter-sized planets (Matsumura 
& Pudritz 2005, 2006). Jets can still be launched from the well-coupled surface layer 
above the dead zóne (e. g. Li 1996; Campbell 2000). 

Although the MRI may be inactive in the bulk of the disc, there are still alternative 
mechanisms of angular momentum transport. In protostellar discs there are probably at 
least two other mechanisms that might contribute to the accretion torque: density waves 
(Róžyczka & Spruit 1993) and the interaction with other planets in the disc (Goodman 
& Rafikov 2001). An additional mechanisms that has been investigated recently is the 
so-called streaming instability that results from the dust trying to move at Keplerian 
speed relative to the surrounding gas that rotates slightly sub-Keplerian due to partial 
pressure support (Youdin & Johansen 2007, Johansen & Youdin 2007). 

Another frequently discussed proposal is the possibility that the turbulence in the 
disc might be driven by a nonlinear hnite amplitúde instability (Chagelishvili et al. 2003). 
While Balbus et al. (1996) and Hawley et al. (1999) have given generál arguments against 
this possibility, Afshordi et al. (2005) and other groups have continued to investigate the 
so-called bypass to turbulence. The basic idea is that successive strong transients can 
maintain a turbulent state in a continuously excited manner. Lesur & Longaretti (2005) 
have recently been able to quantify more precisely the critical Reynolds number required 
for instability. They have also highlighted the importance of pressure fluctuations 
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that demonstrate that a generál argument by Balbus et al. (1996) against nonlinear 
hydrodynamic instabilities is insumcient. 

Yet another possibility is to invoke convectively stable vertical density stratification 
which may give rise to the so-called strato-rotational instability (Dubrulle et al. 2005). 
This instability was recently discovered by Molemaker et al. (2001) and the linear 
stability regime was analyzed by Shalybkov & Rúdiger (2005). However, the presence 
of no-slip radial boundary conditions that are relevant to experiments are vital for the 
mechanism to work. Indeed, the instability vanishes for an unbounded regime, making 
it irrelevant for accretion discs (Umurhan 2006). 

4. Dust dynamics in disc turbulence 

An important aspect of protostellar discs is their capability to produce planets. This 
requires the local accumulation of dust to produce larger conglomerates that can 
condense into rocks and boulders that are big enough to be decoupled from their ambient 
gas flow, and that can combine to grow under the inŕluence of self-gravity. 

The dynamics of dust particles depends critically on their size. Small particles are 
essentially advected by the gas flow, whereas larger ones do not significantly interact 
viscously with the gas and their dynamics is essentially governed by the gravitational 
field leading to Keplerian orbital motion and to settling toward the midplane. If the dust 
grains are small, turbulence will stir up the dust, preventing it from settling. Particles 
of rádius a p can be considered big if its so-called stopping time, r s , is long compared 
with the orbital time. The stopping time is essentially the time it takes for a particle 
to decelerate by a factor of l/e. The precise expression for r s depends critically on the 
ratio of the rádius of the particle, a p , to the mean free path of the gas, t. If a p /i 3> 1, as 
is the case for dust particles in the Earth's atmosphere, one can use the Stokes formula 
for the drag force, 

F d = — QTcpua p u p (Stokes drag), (7) 

where u p is the velocity of the particle. Note that the drag force is linear in u p , so we 
can write F D as iti p u p /t s , where m is the mass of the particle. To get the stopping time, 
this drag force must be balanced by the mass of the particle times its acceleration, i. e. 

m— f = -F D = -m-^. (8) 

dí r s 

We see that the particle speed declines indeed exponentially and that the stopping time 
is given by r s = mu p /Fn, where u p = \u p \. 

In protostellar discs the mean free path is very long, so a p jl <C 1, and one is in the 
so-called Epstein regime were (Seinfeld 1986) 

F Y) Ta —pa^c s u p (Epstein drag). (9) 

This basically means that the mean free path in the kinetic gas formula for the viscosity, 
v = |£c s , is replaced by a p , so the effective viscosity is given by ~ a p c s . 
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Finally, when the velocity of the particle is so large that the Reynolds number with 
respect to the particle rádius, a p u p /u, exceeds a certain critical value, the flow around 
the body is no longer laminar, but turbulent, and so the drag force begins to depend 
quadratically on the velocity of the particle, i. e. 

F r> ~ —pcťpUpU p (turbulent drag). (10) 

In the turbulent regime, the viscosity of the Stokes formula becomes essentially replaced 
by a p u p . 

Let us now estimate the terminál descent speed of a particle in the presence of 
vertical gravity, g = íl 2 z, and let us consider a height of one pressure scale height. Note 
also that in thin discs QH/ y/2 = c s is the sound speed, but in order to avoid additional 
factors of y/2 we prefer to use here H , so that ÍIH = c s . In the following we omit 
the tilde. The terminál descent speed is obtained by balancing the drag force with the 

gravitational force, mg. In the Epstein regime we have therefore 
mg 

u p = — - — (terminál speed). (11) 

P%c s 

Writing this in terms of the Mach number, u p /c s , and expressing the mass in terms of 
rádius and density of the solid, p p ľ^> p, via m ~ p p ap, we have 



Up _ Ppflp g _ p p a p z _ p p a p t z 



(12) 



c s pa 2 p c 2 p H H p t H H 

The mean free path can be expressed in terms of viscosity by using a formula from 
kinetic gas theory, v = |c s £. Thus, £/H = 3u/(c s H) « Re -1 . Using Re = 10 9 for the 
Reynolds number in protostellar discs at 1AU (Brandenburg Sz Subramanian 2005), 
together with z ~ H, and p p /p ~ 10 10 (see Hodgson & Brandenburg 1998), we have 

£ « 10 x °f. (13) 

The dependence u p /c s versus a p /£ is shown in Fig. [2] both for laminar flows in the 
Epstein regime and for turbulent flows where the product of u p /c s and a p /£ exceeds 
unity. The stopping time expressed in non-dimensional form is given by 

mup ppdp 

Qr s = V— = --. (14) 
Using again p p /p ~ 10 10 we find ílr s = 1 for meter-sized particles. 

4-1. Particle trapping in vortices 

In the following we discuss why in the size range where Qt s = 1 particles can be trapped 
by vortices. This possibility was first suggested by Barge & Sommeria (1995) and 
Tanga et al. (1996) who proposed that dust particles could be trapped within turbulent 
anticyclonic eddies. This mechanism can be understood by taking the divergence of the 
evolution equations for the particle velocity, 

du n — . 2 _ x 



dt 



g - 2fž x Up + 3fi Xp - t s 1 (u p - u), (15) 
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Figúre 2. Regimes of applicability of Stokes, Epstein, and turbulent drag formulac 
as a funetion of particle rádius, a p , and particle speed, u p . The dotted lines gives the 
terminál descent speed of particles as a funetion of a p . Note that in the turbulent 
regime the slope of the dotted eurve is 1/2, whilc in the Epstein regime it is 1. For 
comparison, l w 10 m for the solar nebula at 1 AU. 



where u is the velocity of the gas. Taking separately the curl and the divergence of this 
equation we obtain two scalar equations, 

^ (20 • Wp ) = -Aíl 2 (V • u p ) - 2QT- 1 • (u> p - w), (16) 

j t (V • u p ) = 2íl-u> p + 2tf - t-\V ■ u p ), (17) 

where u> = V x u and u p = V x u p are the vorticity of gas and particles, respectively, 
and we have assumed that the gas is solenoidal, i. e. V • u, and that only horizontál 
velocity components enter. Eliminating 2fl-uj p from these two equations yields (Hodgson 
& Brandenburg 1998) 

[(2Qr s ) 2 + (1 + r s d t ) 2 } (V • u p ) = 2Slr s ■ (w + íl) . (18) 

This equation shows that there will be a trend toward a negative divergence (i. e. a 
positive convergence) of the particle velocity if 

2ÍI ■ (w + íl) < 0. (19) 

Thus, for particle aceumulation, i. e. V • u p < 0, not only does the vorticity have to 
be anticyclonic, but also the vorticity has to be sufnciently anticyclonic so that the 
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Figúre 3. Gas and dust in the mid-plane after one orbít. Gas density and velocity 
field for the three values of ŕŽT s are indistinguishable, so only fŽT s = 1 is shown. The 
other three plots show dust density and dust velocity field after one orbit for fžr s = 0.1, 
ížt s = 1 and S!t s = 10 respectively. The dust velocity field of fžr s = 0.1 is very similar 
to that of the gas, due to the short stopping time. For Qt s = 1 there is a strong 
convergence towards the interior of the vortex, whereas for ížt s = 10 only a slight 
density inerease in a narrow región that extends from the vortex along the shear is 
seen. Courtesy of Johansen et al. (2004). 

condition above is satisfied. Looking at histograms of vorticity, it is clear that this 
condition is only satisfied in the very tail of the distribution of the axial component 
of the vorticity (see Fig. 2 of Hodgson & Brandenburg 1998). The trapping of gas by 
anticyclonic vortices has been studied by a number of people (Barge & Sommeria 1995, 
Tanga et al. 1996, Klahr & Henning 1997, Chavanis 2000, Johansen et al. 2004). In Fig. [3] 
we reproduce the result of a simulation by Johansen et al. (2004), where density and 
velocity of the gas are compared with density and velocity of the dust for an anticyclonic 
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Figúre 4. Simulations of Johansen et al. (2007) showing the total column density (gas 
plus all particle sizes) in the horizontál piane. The insets show the column density in 
logarithmic scale centered around the most massive cluster in the simulation. Time 
is given in terms of the number of orbits after turning on self-gravity. Courtesy of 
Johansen et al. (2007). 



vortex with three values of fžr s . 

4-2. Planetesimal formatíon with self-gravity 

A number of new simulations have emerged in recent years. A major step was made in 
a paper by Johansen et al. (2007) who combined the dust dynamics with self-gravity in 
the shearing box approximation (Fig. 4). One of the remarkable results they find is a 
rapid formation of Ceres-sized bodies from boulders. Even though the mass of what one 
might call protoplanet is growing, this body is also shedding mass during encounters 
with ambient materiál as it flows by One might speculate that what is missing is 
the effect of radiative cooling of the protoplanet. This would allow the newly acereted 
materiál to lose entropy, become denser, and hence fall deeper into its potential well. 

5. Outflows from protostellar discs 

Finally we diseuss the phenomenon of outflows and collimated jets from protostellar 
discs. A common approach to modelling jets is by treating the disc as a boundary 
condition where a poloidal magnetic field is wound up by Keplerian rotation in the disc. 
A particularly useful setup has been diseussed and studied in the papers by Ouyed & 
Pudritz (1997a,b) and Fendt & Elstner (1999). 

By modelling the disc as a boundary condition, it is impossible to account for 
the generation and evolution of magnetic fields that are responsible for the centrifugal 
acceleration of outflows via the Blandford and Payne (1982) mechanism. This was the 
reason why von Rekowski et al. (2003) proposed a mean field dynamo model that allows 
for outflows. We briefly summarize some of the main findings. 
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Figúre 5. General structure of the outflows typically obtained in the model of von 
Rekowski et al. (2003), where the cool, dense disc emits a thermally driven wind and a 
magneto-centrifugally driven outflow in a conical shell. [Adapted from von Rekowski 
et al. (2003).] 



5.1. Models wíthout magnetízed centrál star 

A mean field dynamo model with piece-wise polytropic hydrodynamics was proposed 
by von Rekowski et al. (2003). Like Ouyed & Pudritz (1997a,b) they štart with an 
equilibrium corona, where they assume constant entropy and hydrostatic equilibrium 
according to c p T(r) = GM*/r, i. e. enthalpy equals the negative gravitational potential. 
In order to make the disc cooler, a geometrical región for the disc is prescribed (see 
Fig. [5]). The entropy contrast between disc and corona is chosen such that the initial 
disc temperature is about 3 x 10 3 K in the bulk of the disc. The low disc temperature 
corresponds to a high disc density of about 10 -10 . . . 10 _9 gcm -3 . For the disc dynamo, 
the most important parameter is the dynamo coefficient ttdyn in the mean field induction 
equation. The dynamo a effect is antisymmetric about the midplane and restricted to 
the disc. 

The upper panel of Fig. [6] illustrates that an outflow develops that has a well- 
pronounced structure. Within a conical shell originating from the inner edge of the disc 
the terminál outflow speed exceeds 500 km/s and temperature and density are lower 
than elsewhere. The inner cone around the axis is the hottest and densest región where 
the stellar wind speed reaches about 150 km/s. The wind that develops from the outer 
parts of the disc has intermediate values of the speed. 

The structured outflow is driven by a combination of different processes. A 
significant amount of angular momentum is transported outwards from the disc into 
the wind along the magnetic field, especially along the concentrated lines within the 
conical shell (see the lower panel of Fig. [6]). The magnetic field geometry is such that 
the angle between the rotation axis and the field lines threading the disc exceeds 30° at 
the disc surface, which is favorable for magneto-centrifugal acceleration (Blandford Sz 
Payne 1982). However, the Alfvén surface is so close to the disc surface at the outer parts 
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Figúre 6. Left: poloidal velocity vectors and poloidal magnetic field lines (white) 
superimposed on a color scale representation of log h. Specific enthalpy h is directly 
proportional to temperature T, and logh — (—2,-1,0,1) corresponds to T w 
(3xl0 3 ,3xl0 4 ,3xl0 5 ,3xl0 6 )K. The black dashed Ľne shows the fast magnetosonic 
surfacc. The disc boundary is shown with a thin black line, the stellar surface is marked 
in red. The dynamo adyn coefíicient is negative in the upper disc half, resulting in 
roughly dipolar magnetic symmetry. Right: color scale representation of the specific 
angular momentum, normalized by the maximum angular momentum in the disc, with 
poloidal magnetic field lines superimposed (white). The black solid line shows the 
Alfvén surface, the white dashed line the sonic surface. Same model as in the left hand 
panel and averaged over same times. [Adapted from von Rckowski & Brandenburg 
(2004).] 

of the disc that acceleration there is mainly due to the gas-pressure gradient. In the 
conical shell, however, the outflow is highly supersonic and yet sub-Alfvénic, with the 
Alfvén rádius a few times larger than the rádius at the footpoint of the field lines at the 
disc surface. The lever arm of about 3 is sufficient for magneto-centrifugal acceleration 
to dominate in the conical shell (cf. Krasnopolsky et al. 1999). 

5.2. Star-disc couplíng 

We now diseuss the interaction of a stellar magnetic field with a circumstellar aceretion 
disc and its magnetic field. This problém was originally studied in connection with 
aceretion discs around neutrón stars (Ghosh & Lamb 1979), but it was later also applied 
to protostellar magnetospheres (Kônigl 1991, Cameron & Campbell 1993, Shu et al. 
1994). Most of the work is based on the assumption that the field in the disc is constantly 
being dragged into the inner parts of the disc from large radii. The underlying idea is 
that a magnetized molecular cloud collapses, in which case the field in the centrál star 
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Figúre 7. Sketch showing the formation of an X-point when the disc field is aligned 
with the dipóle (on the left) and the formation of current sheets with no X-point if 
they are anti-aligned (on the right). The two current sheets are shown as thick lines. 
In the present paper, the second of the two configurations emerges in all our models, 
i. e. with current sheets and no X-point. [Adapted from von Rekowski & Brandenburg 
(2004).] 

and that in the disc are aligned. This scenario was studied numerically by Hirose et 
al. (1997) and Miller & Stone (1997). In the configurations they considered, there is 
an X-point in the equatorial piane (see left hand panel of Fig. [7j), which can lead to a 
strong funnel flow. 

However, the opposite relative orientation is in principle also possible and has been 
explored by Lovelace et al. (1995), where the magnetic field of the star has been flípped 
and is now anti-parallel with the field in the disc, so that the field in the equatorial 
piane points in the same direction and has no X-point. However, current sheets develop 
above and below the disc piane (see right hand panel of Fig. [7]). This configuration is 
also referred to as the X-wind model. Ironically, this is the field configuration without 
an X-point. 

Simulations of such a field configuration by Hayashi et al. (1996) confirm the idea 
by Lovelace et al. (1995) that closed magnetic loops connecting the star with disc are 
twisted by differential rotation between star and disc, and then inflate to form open 
stellar and disc field lines (see also Bardou 1999, Agapitou & Papaloizou 2000, Uzdensky 
et al. 2002). Goodson et al. (1997, 1999) and Goodson & Winglee (1999) find that for 
sufficiently low resistivity, an accretion process develops that is unsteady and proceeds 
in an oscillatory fashion. Similar results have also been obtained by von Rekowski & 
Brandenburg (2004). In their case the resulting field geometry is always the one in the 
second panel of Fig. [3, i. e. the one with current sheets and no X-point. It is seen that 
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the inflating magnetosphere expands to larger radii where matter can be loaded onto 
the field lines and be ejected as stellar and disc winds. Furthermore, reconnection of 
magnetic field lines allows matter to flow along tliem and accrete onto tlie protostar, in 
tlie form of a funnel flow (Romanova et al. 2002). 

6. Conclusions 

Protostellar discs are generally believed to be turbulent, although there is substantial 
uncertainty regarding tlie relative importance of possible mechanisms and tlie necessary 
physics involved. The turbulence in protostellar discs may not be as vigorous as in fully 
ionized discs around neutrón stars and black holes, but its presence is critical for driving 
accretion toward tlie centrál star, driving mass concentrations within tlie disc above the 
critical value for gravitational collapse to form planets. Turbulence is also important 
for heating a corona above the disc from where outflows can be driven. 

Several important aspects have been omitted in this review. The inclusion of 
radiative cooling has already been mentioned as one of the important ingredients that 
allow matter to settle deeper in the potential well and hence to prevent newly accreted 
materiál from being stripped away from the forming protoplanet. Related to this is 
the question how effective the acceleration of outflows is when one allows for radiative 
cooling of the jet. Another aspect concerns the generation of large scale magnetic fields 
both in the protostar and the disc. Much attention has been paid to the importance of 
shedding small scale magnetic helicity from the dynamo, because otherwise the dynamo 
will suffocate from excess small scale magnetic helicity that quenches the production of 
large scale magnetic fields by kinetic helicity and/or shear. Advanced global simulations 
are likely to shed light on these and related questions. 

References 

[I] Afshordi, N., Mukhopadhyay, B., Narayan, R. 2005, ApJ, 629, 373 
[2] Agapitou, V., & Papaloizou, J. C. B. 2000, MNRAS, 317, 273 

[3] Armitage, P. J. 1998, ApJ, 501, L189 

[4] Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 

[5] Balbus, S. A. & Hawley, J. F. 1998, Rev. Mod. Phys., 70, 1 

[6] Balbus, S. A., Hawley, J. F., & Stone, J. M. 1996, ApJ, 467, 76 

[7] Bardou, A. 1999, MNRAS, 306, 669 

[8] Bardou, A., Rekowski, B. v., Dobler, W., Brandenburg, A., et al. 2001, A&A, 370, 635 
[9] Barge, P., & Sommeria, J., 1995, A&A, 295, Ll 

[10] Blackman, E. G., Penna, R. F., & Varniérc, P. 2008, New Astron., 13, 244 

[II] Blandford, R. D., & Payne, D. R. 1982, MNRAS, 199, 883 

[12] Brandenburg, A. 1998, in Thcory of Black Hole Accretion Discs, ed. M. A. Abramowicz, G. 

Bjornsson & J. E. Pringle (Cambridge University Press), 61 
[13] Brandenburg, A., & Subramanian, K. 2005, Phys. Rep., 417, 1 
[14] Brandenburg, A., Nordlund, Á, Stcin, R. F., et al. 1995, ApJ, 446, 741 
[15] Brandenburg, A., Nordlund, Á, Stein, R. F., et al. 1996, ApJ, 458, L45 
[16] Cameron, A. C, & Campbell, C. G. 1993, A&A, 274, 309 



Turbulent protostellar discs 



[17] Campbell, C. G. 2000, MNRAS, 317, 501 

[18] Chagelishvili, G. D., Zahn, J.-R, Tevzadze, A. G., & Lominadzc, J. G. 2003, A&A, 402, 401 

[19] Chavanis, P. H., 2000, A&A, 356, 1089 

[20] Dubrullc, B., Márie, L., Normand, et al. 2005, A&A, 429, 1 

[21] Fendt, C, & Elstner, D. 1999, A&A, 349, L61 

[22] Fleming, T., & Stone, J. M. 2003, ApJ, 585, 908 

[23] Fromang, S., Tcrqucm, C, & Balbus, S. A. 2002, MNRAS, 329, 18 

[24] Fromang, S., & Papaloizou, J. 2007, A&A, 476, 1113 

[25] Fromang, S., Papaloizou, J., Lesur, G., & Heinemann, T. 2007, A&A, 476, 1123 

[26] Fromang, S., & Nclson, R. P. 2006, A&A, 457, 343 

[27] Gammie, C. F. 1996, ApJ, 457, 355 

[28] Ghosh, P. & Lamb, F. K. 1979, ApJ, 232, 259 

[29] Glassgold, A. E., Najita, J., & Igea, J. 1997, ApJ, 480, 344 

[30] Goodman, J., & Rafikov, R. R. 2001, ApJ, 552, 793 

[31] Goodson, A. P., & Winglec, R. M. 1999, ApJ, 524, 159 

[32] Goodson, A. P., Bohm, K.-H., & Winglee, R. M. 1999, ApJ, 524, 142 

[33] Goodson, A. P., Winglee, R. M., & Bohm, K.-H. 1997, ApJ, 489, 199 

[34] Hawley, J. F. 2000, ApJ, 528, 462 

[35] Hawley, J. F., Gammie, C. F., & Balbus, S. A. 1995, ApJ, 440, 742 

[36] Hawley, J. F., Gammie, C. F., & Balbus, S. A. 1996, ApJ, 464, 690 

[37] Hawley, J. F., Balbus, S. A., & Winters, W. F. 1999, ApJ, 518, 394 

[38] Hayashi, M. R., Shibata, K., & Matsumoto, R. 1996, ApJ, 468, L37 

[39] Hirose, S., Uchida, Y., Shibata, K., et al. 1997, Publ. Astron. Soc. Japan, 49, 193 

[40] Hodgson, L. S. & Brandcnburg, A. 1998, A&A, 330, 1169 

[41] Igea, J., & Glassgold, A. E. 1999, ApJ, 518, 848 

[42] Inutsuka, S.-L, & Sano, T. 2005, ApJ, 628, L155 

[43] Johanscn, A., Andcrscn, A. C, & Brandenburg, A. 2004, A&A, 417, 361 

[44] Johanscn, A., & Youdin, A. 2007, ApJ, 662, 627 

[45] Johanscn, A., Oishi, J. S., Mac Low, M. M., et al. 2007, Náture, 448, 1022 

[46] Klahr, H. H. & Hcnning, T. 1997, Icarus, 128, 213 

[47] Konigl, A. 1991, ApJ, 370, L39 

[48] Krasnopolsky, R., Li Z.-Y., & Blandford R. D. 1999, ApJ, 526, 631 

[49] Lesur, G., & Longaretti, P.-Y. 2005, A&A, 444, 25 

[50] Li, Z.-Y. 1996, ApJ, 465, 855 

[51] Lovclace, R. V. E., Romanova, M. M., & Bisnovatyi-Kogan, G. S. 1995, ApJ, 275, 244 

[52] Lyra, W., Johanscn, A., Klahr, H., & Piskunov, N. 2008, A&A, 479, 883 

[53] Matsumura, S., & Pudritz, R. E. 2005, ApJ, 618, L137 

[54] Matsumura, S., & Pudritz, R. E. 2006, MNRAS, 365, 572 

[55] Miller, K. A., & Stone, J. M. 1997, ApJ, 489, 890 

[56] Molemaker, M. J., McWilliams, J. C, & Yavnch, I. 2001, Phys. Rev. Lctt., 86, 5270 

[57] Nclson, R. P., & Papaloizou, J. C. B. 2003, MNRAS, 339, 993 

[58] Oishi, J. S. Mac Low, M. M., Menou, K. 2007, ApJ, 670, 805 

[59] Ouyed, R., & Pudritz, R. E. 1997a, ApJ, 482, 712 

[60] Ouyed, R., & Pudritz, R. E. 1997b, ApJ, 484, 794 

[61] Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., et al. 2004, ApJ, 616, L151 

[62] Róžyczka, M., & Spruit, H. C. 1993, ApJ, 417, 677 

[63] Sano, T., & Stone, J. M. 2002, ApJ, 570, 314 

[64] Sano T., Miyama S. M., Umebayashi T., and Nakano T. 2000, ApJ, 543, 486 

[65] Seinfeld, J. H. 1986, Atmosphcric Chcmistry and Physics of Air Pollution (Wiley, New York) 

[66] Semenov, D., Wiebe, D., & Hcnning, T. 2004, A&A, 417, 93 

[67] Shakura, N. L, & Sunyaev, R. A. 1973, A&A, 24, 337 



Turbulent protostellar discs 



16 



[68] Shalybkov, D. & Riidiger, G. 2005, A&A, 438, 411 

[69] Shu, F., Najita, J., Ostriker, E., et al. 1994, ApJ, 429, 781 

[70] Steinacker, A., & Papaloizou, J. C. B. 2002, ApJ, 571, 413 

[71] Stone, J. M., Hawley, J. F., Gammie, C. F., et al. 1996, ApJ, 463, 656 

[72] Tanga, P., Babiano, A., Dubrulle, B., & Provenzale, A. 1996, Icarus, 121, 158 

[73] Torkelsson, U., Brandenburg, A., Nordlund, Ä., et al. 1996, Astrophys. Letter & Comm., 34, 383 

[74] Umurhan, O. M. 2006, MNRAS, 365, 85 

[75] Uzdcnsky, D. A., Konigl, A., & Litwin, C. 2002, ApJ, 565, 1191 

[76] von Rekowski, B., & Brandenburg, A. 2004, A&A, 420, 17 

[77] von Rekowski, B., Brandenburg, A., Dobler, W., et al. 2003, A&A, 398, 825 

[78] Youdin, A., & Johansen, A. 2007, ApJ, 662, 613 



